The Milky Way's stellar halo - lumpy or triaxial? 
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' Abstract. We present minimum chi-squared fits of power law and Hernquist density profiles 

to F-turnoff stars in eight 2.5° wide stripes of SDSS data: five in the North Galactic Cap 
and three in the South Galactic cap. Portions of the stellar Galactic halo that are known to 
contain large streams of tidal debris or other lumpy structure, or that may include significant 
contamination from the thick disk, are avoided. The data strongly favor a model that is not 
symmetric about the Galaxy's axis of rotation. If included as a free parameter, the best fit to 
the center of the spheroid is surprisingly « 3kpc from the Galactic center in the direction of the 
Sun's motion. The model fits favor a low value of the density of halo stars at the solar position. 
■ The alternative to a non-axisymmetric stellar distribution is that our fits are contaminated by 

previously unidentified lumpy substructure. 
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1. Introduction 

The stars in the Milky Way are divided into: an old, metal-rich bulge population near the 
Galactic center; a thin disk population in the plane of the Milky Way that includes a range of 
stellar ages and metallicities but is unique in including a young, metal-rich population; a thick 
^ ■ disk population of intermediate age and metallicity (fraction of elements heavier than helium); 

and an old, metal-poor spheroid (also called the stellar halo) population in which the disks are 
embedded. Each of the populations has a distinct kinematic signature. The spheroid is generally 
described as a somewhat spherical population with a density that falls off as r _Q: , a = 3.5, where 
r is the distance from the center of the Galaxy. A squashing parameter 0.5 < q < 1.0, in the 
direction perpendicular to the Galactic plane is also fit to the density distribution of the spheroid 
by redefining the distance from the Galactic center as: r = \/X 2 + Y 2 + (Z/q) 2 . 

The Sloan Digital Sky Survey (SDSS) scans the sky in 2°. 5 wide stripes that follow great 
circles on the sky. The photometric images are taken in five optical filters: ugriz. Data that has 
been corrected for interstellar reddening using IRAS dust maps [I] is designated by a subscripted 
zero. Three public data releases [21 El HI are available though an online catalog archive server. 
Though the survey geometry and technical design [SI El El El El QUI CD] optimized extragalactic 
research, studies of the Galactic stars in the photometric catalog have significantly advanced 
our knowledge of the stars in the halo of our own Galaxy. 

In particular, it has been shown 12 that the distribution of faint F-turnoff stars (19 < 
g < 21.5) that we generally associate with the spheroid population [121 EI are n °t distributed 
symmetrically about the Galaxy center. Even if regions surrounding known or suspected tidal 
streams and debris are eliminated from consideration, the stellar distribution is not symmetric 
about a plane (perpendicular to the Galactic plane) through Galactic longitudes / = 0°, / = 180°. 
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Figure 1. Location of the DR3 data used to fit the spheroid in Galactic coordinates. The grey 
scale shows spheroid stars selected from DR3 in the manner described in 12 J. The black lines 
show the outlines of eight stripes used in this paper to fit the spheroid models. The SDSS stripes 
numbered are sequentially in a clockwise direction around the coordinate (l,b) = (209°,— 7°). 
This figure shows stripes 10 (bottom) and 15 (top) in the range 0° < / < 29°, 209° < / < 360°; 
and stripes 24, 32, 37, 76, 82, 86, top to bottom, in the range 29° < / < 209°. 



Since there are no values of the axisymmetric power law model parameters (a, q) which describe 
an asymmetric distribution, the standard model must be modified to include asymmetric stellar 
populations, large scale lumpiness, or an asymmetric smooth spheroid component. 

In this paper, we make the assumption that there exists a smooth component to the Galactic 
halo, and present several candidate models. We then use a minimum x 2 technique to find the 
best fit parameters for each model, fitting to eight 2°. 5 wide stripes of data in the SDSS DR3. 
These stripes are chosen to sample the widest range of positions available in the Galaxy from 
this dataset. Figure [l] shows the positions of these data stripes in Galactic coordinates. 

2. Spheroid Models and Parameters 

The simplest extension to the oblate power law spheroid model that allows for asymmetry about 
the Galactic center is a triaxial power law spheroid. In addition to flattening the spheroid by 
a factor of q in the Z direction, we flatten by a factor of p in the Y direction. This density 
function with ellipsoidal contours can in general be shifted in three coordinates and rotated by 
three angles with respect to the Galactic coordinates (X, Y, Z), where X goes from the solar 
position through the center of the Galaxy, Y is in the direction of the local standard of rest 
(roughly the direction the Sun is going), and we have a right-handed coordinate system centered 
on the center of the Milky Way. The principle axes of the density function are: 

X r = x cos 9 cos (p + y sin 9 cos (p + z sin (p 

Y f = x(— sin 9 cos £ — cos 9 sin (p sin £) + 2/(cos#cos£ — sin 9 sin <p sin £) + z cos </> sin £ 
Z' = x(sin^sin^ — cos 9 sin (p cos £) + y(— cos 9 sin £ — sin 9 sin (p cos £) + zcos^cos^, 



where x = X — dX, y = Y — dY, and z = Z — dZ. Initially, we expected that dX dY = dZ = 
= £ = 0, but the addition of these parameters to the models significantly improved the fits. 

We also found that power laws tended to predict too high a star density near the center of 
the Galaxy, compared to our data. We achieved a better fit to the data with a Hernquist profile: 

P = r{r + X ore )3 > r = y/x* + ( y Vp) 2 + 
Although we did explore the possibility of fitting a function of the form: 



H r<*(r + R core y 

which would include both the power law model (5 = 0) and the Hernquist model (a = 1, S = 3), 
this fitting was abandoned because we are not very sensitive to the combination of a and S. In 
general, models with a + 5 w 4 worked well. 



3. F Turnoff Star data 

The data includes all stellar objects in the 8 SDSS stripes with: (r — z)o < 0.8 (to get rid of 
M stars at faint magnitudes), 0.4 < (u — g)o < 3.8 (to eliminate QSOs), 0.0 < (g — r)o < 0.3 
(to select stars with turnoff colors), and go < 24. The designation as "star" in the database 
removes issues of duplication of objects due to overlaps and deblending. Objects are further 
rejected if flags indicate that they are saturated, too near the edge of an image to give good 
photometry, or classified as a "bright" object. This resulted in selection of (234,306; 190,221; 
159,632; 132,637; 123,359; 115,416; 165,620; 150,562) stars from stripes (10, 15, 24, 32, 37, 76, 
82, 86), respectively. These stripes are shown in figured 

In order to compare the stellar density to the model density for a given set of parameters, the 
stellar data was sorted into bins that were 0.1 mag. wide in go and 1° in angle along the length 
of the stripe. For a given model, the expected number of stars in each bin was calculated as: 

N = p(r)dV = p(r)R 3 , 
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where R is the distance from the Sun, dm is the magnitude bin size, is the solid angle of each 
bin (2.5 sq. deg.), and p is the density of our sample of F turnoff stars. 

Since our models included only stars in a smooth spheroid population, we selected data only 
from parts of sky and magnitude ranges we thought were most likely to contain exclusively 
stars from this population. Figures El El El 13 El El and El show the outlines of the portions 
of the data that seemed the most likely to contain purely spheroid stars. Bright magnitudes 
were avoided, particularly near the Galactic plane, since they could be contaminated with thick 
disk stars (as determined by Besancon galaxy models [El)- Stars near any previously identified 
spheroid over density were also avoided. 

The models were then modified to reflect the uncertainty in the absolute magnitudes of the 
turnoff stars, which leads to an uncertainty in the distance to each star. We analyzed the stars 
in the three globular clusters Pal 5 (stripe 10), M15 (stripe 76), and M2 (stripe 82) to determine 
the turnoff star counts as a function of apparent magnitude. A histogram of star counts vs. go 
magnitude was created by subtracting the star counts in a neighboring piece of sky from the 
star counts in the direction of the globular cluster. A function of the form: 

f (an) - A p -(90-9 P eak)/C2a 2 ) f O- = 0.55, go < g P eak 

J[90) I (T = 0.75 , 9 >g P eak 



was a reasonable fit to the data. This is an approximate fit; each of the clusters had a slightly 
different magnitude profile. 



So that the model took into account this range of absolute magnitudes, we convolved the 
model with this magnitude distribution, adjusting the amplitude A so that the sum of the bins 
in the convolution kernel was 1.0. This way, the overall normalization of the model was not 
affected. 

4. Minimum x 2 fits 

Once a set of data was extracted and a set of models was defined, we fit the models to the data 
by adjusting the model parameters to produce the minimum x 2 - We considered five models: a 
triaxial Hernquist model in which the center was not fixed to the center of the Galaxy, a triaxial 
Hernquist model in which the center was forced to the center of the Galaxy, a triaxial power law 
model in which the center was allowed to float, a triaxial power law model in which the center 
was forced to the center of the Galaxy, and a traditional axisymmetric power law model. The 
models contained between 4 and 10 free parameters, depending on the model. 

For each model, we chose a starting set of parameters, excluding the normalization po 5 an d 
generated the expected star counts for each bin in the star count data. We then normalized the 
model so that the sum of all of the star counts equaled the sum of all of the star counts in the 
data set (including all stripes) that we were fitting. This had the advantage of reducing by one 
the number of parameters that were allowed to freely vary, since this one could be computed in 
a straightforward manner. The value of x 2 f° r this model and this parameter set is then: 

2= E?=b n (d l -m l ) 2 /m l 

■^bin ^param 2 

where Nu n is the number of bins used in all stripes, d{ represents each bin of the data, rrii 
represents the expected number of counts in that bin in the model, and N param is the number 
of parameters fit in the model. We have assumed that the statistical errors, cr^, in the number 
of stars in the bins are Poisson so that of = m^. 

We changed the parameters in each model by hand until we achieved the minimum x 2 f° r 
each model. The final parameters were checked to guarantee that the fit is worse if any of the 
parameters are changed by ± the tolerance (given in column 2 of table CJ for that parameter. If 
there are 10 adjustable parameters, then 3 10 parameter sets were checked to ensure a minimum. 

The best fit model parameters are given in tabled Parameters include: i?o, distance from 
the Sun to the Galactic center; (p 5 q), the ratio of the scale length in the (Y f ,Z f ) directions 
compared to the major axis in the X' direction; R CO rei the Hernquist core radius; Mj, the peak 
of the absolute magnitude distribution of our turnoff stars; and po 5 the normalization of the 
Hernquist or power law model. (a,S) are (1,3) for a Hernquist profile and (a, 0) for a power 
law profile. To get from Galactic XYZ coordinates to the principle axes of the triaxial models, 
we first shift the center by (dX,dY,dZ), then rotate by 9 around the Z-axis, then rotate by 
(j) around the new F-axis, and then by £ around the new X-axis. N so i ar is the density of the 
selected turnoff stars in the solar neighborhood and x 2 5 & measure of the goodness of the model 
fit. 

Numbers in bold were set to a fixed value and were not allowed to vary. Parameters were 
fixed either because the definition of the model required it (for example, power laws all have 
8 — 0), or because the parameter was correlated with a parameter that was varied. For instance, 
in the full Hernquist model we set the scale of the model with R CO re, though we could instead 
have varied Mf. Ro and dX are redundant. 

Figures El El El El El HI El and El show the star density in the data and the residual star 
densities in each stripe of data after the model has been subtracted. Figure El additionally shows 
the star density of the model for comparison. 

A crude estimate for the mass density of stars in the solar neighborhood is calculated from 
SDSS measurements of the Pal 5 globular cluster. There are 588 stars above background with 



Table 1. Best Fit parameters for five spheroid models described in Section 2. Numbers in bold 
were fixed; all others were varied to minimize x 2 fit between the model and the data The units 
of po are kpc 4 for Hernquist models, and kpc a for power law models. 
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0.0 < g -r < 0.3 in Pal 5. The mass of Pal 5 is 4.5 to 6.0 x 1O 3 M [TH]. Assuming that Pal 5 is 
a good approximation to the stellar population of the spheroid (a big assumption), an estimate 
of the number of solar masses per sample star is lOM^/star. Therefore, the local mass density 
of the stellar spheroid is 11,000 M /kpc 3 for the full Hernquist model, and 15,000 M /kpc 3 
for the best fit standard power law. These numbers are 20-25 times lower than the average 
density from Allen's Astrophysical Quantities of 2.6 x lO 5 M /kpc 3 . Since Pal 5 has a bluer 
turnoff than the spheroid, it is likely that using this cluster to estimate the mass of the spheroid 
component leads to an underestimate, though it seems unlikely that it could account for the 
whole difference. Future work will provide a better estimate of the spheroid normalization. 

5. Conclusions 

The standard stellar density profile used to fit the Galactic spheroid, p = po/r a , is an extremely 
poor fit to the density of halo turnoff stars when a large area of the sky is sampled, even using 
the unusually large best fit distance between the Sun and the Galactic center. The flattening of 
q = 0.63 and the power law slope of a = 3.1 are within the range of previously measured values. 

Triaxial models are best fit with a major axis 50 — 70° from the line of sight from the Sun 
to the center of the Galaxy and tilted 4° — 6° from the Galactic plane. The minor axis, ^ 13° 
from the Z-direction, is about 65% of the scale length of the major axis. The intermediate axis 
is about 75% of the scale length of the major axis, and tilted at about 12° from the Galactic 
plane. These models are clearly favored over the axisymmetric standard model. 

The Hernquist profile is a better fit to the data than a power law profile, but give roughly 
similar values for the positional and flattening parameters. All of the models suggest a lower 
spheroid stellar density near the Sun compared with previous estimates. 

If the center of the spheroid profile is allowed to vary, then it moves a surprising 3.0 - 3.5 




Figure 2. Stripe 37 model, data, and subtraction. The radial coordinate is the reddening 
corrected go magnitude; if we assume the turnoff stars all have M go = 4.2, then the radial 
coordinate is logarithmically proportional to distance from the Sun. The approximate distances 
in kpc are given on the bottom half of the model and data plots. The model and the data are 
shown with the same grey scale. The two outlined regions on the data and the subtractions 
show the portions of the data on this stripe that were used to fit the model. We avoided the 
faintest stars because they might not be complete, the brightest stars because they could be 
contaminated by thick disk stars, and the portion of the data at the top of the diagram that 
includes a known overdensity in the Galactic plane ^1]. The scale for the model and data is 
such that zero counts is white; in the subtraction zero counts is grey so that negative deviations 
(white) and positive deviations (black) are both discernable. 

kpc from the nominal center of the disk in the direction of the Sun's motion. Since this is 
perpendicular to the line of sight from the Sun to the Galactic center, it is not possible to 
reduce this offset by simply changing the scale lengths or assumed absolute magnitudes of the 
sample stars. It is the stars in the southern stripes near the Galactic center that tend to pull 
the spheroid density off center. We are still investigating whether there could be unidentified 
substructure or data quality issues at the root of this surprising result. Clearly, more data in 
the southern hemisphere would be useful in understanding the large-scale distribution of stars 
in the halo. 

A large fraction of the Galactic halo contains known substructure that must be avoided to 
fit the smooth component. It is unknown whether the results provided here are contaminated 
by unidentified substructure. 

Acknowledgments 

H. J. N. acknowledges funding from Research Corporation and the National Science Foundation 
(AST 03-07571). Ken Freeman and Markus Samland provided stimulating discussion of the of 
the stellar content of the Milky Way. 

Funding for the creation and distribution of the SDSS Archive has been provided by the 
Alfred P. Sloan Foundation, the Participating Institutions, the National Aeronautics and Space 
Administration, the National Science Foundation, the U.S. Department of Energy, the Japanese 



[223,20] [223,20] 




Figure 3. Stripe 10 data and subtraction. Very little of this stripe could be used for comparison 
due to the tidal stream in the plane near the top of the stripe, the Sagittarius stream on the 
bottom left, and an apparent overdensity of unknown origin on the middle left. The outlined 
region on the top of the stripe was included because some of the model fits tended to oversubtract 
this region. 
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Figure 4. Stripe 15 data and subtraction. Prominent again in this figure are the Sagittarius 
dwarf tidal stream (left) and a tidal stream in the Galactic plane (top). The excess counts in the 
subtraction near the bottom of the wedge, near the Galactic center, could include stars leaking 
into the sample from the thick disk. 




Figure 5. Stripe 24 data and subtraction. This stripe has fewer large tidal streams. The stream 
in the Galactic plane is still present at the top of the diagram, and there is a small overdensity 
of unknown origin about 30° away in the counterclockwise direction. Again, brighter stars near 
the Galactic center were avoided in the fitting procedure due to possible contamination from 
the thick disk. 




Figure 6. Stripe 32 data and subtraction. This is a very clean stripe, with a nice subtraction. 
The stream in the plane of the Milky Way is present at the top of the figure. 
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Figure 7. Stripe 76 data and subtraction. This is a shorter section of data that is free from 
tidal debris, but contains a section of poor data and a globular cluster (M15) that appear as 
dark radial features. 




Figure 8. Stripe 82 data and subtraction. The trailing tail of the Sagittarius tidal stream is 
prominent here. Also, there is a section of poor data and a globular cluster (M2) that show up 
as radial features most easily seen in the subtraction. The poor subtraction near the Galactic 
center (bottom) at bright magnitudes could be due to thick disk stars leaking into the sample. 




Figure 9. Stripe 86 data and subtraction. The Sagittarius stream trailing tidal tail is also 
prominent in this stripe. The data near the Galactic center shows a radial pattern that is 
indicative of poorer photometric data. We experimented with leaving this stripe in and taking it 
out. In the end we decided to keep the data since the fits to this stripe were fairly consistent with 
the other southern stripes, and keeping it improved the overall appearance of the subtractions. 
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